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JAN MANSCHOT 

Abstract. Generating functions of BPS invariants for JV = 4 U(r) gauge theory on a Hirze- 
bruch surface with r < 3 are computed. The BPS invariants provide the Betti numbers of 
moduli spaces of semi-stable sheaves. The generating functions for r = 2 are expressed in 
terms of higher level Appell functions for a certain polarization of the surface. The level 
corresponds to the self-intersection of the base curve of the Hirzebruch surface. The non- 
holomorphic functions are determined, which added to the holomorphic generating functions 
provide functions which transform as a modular form. 



1. Introduction 

The study of supersymmetric spectra of field theories and supergravities is a major subject 
in theoretical physics and also mathematics. The BPS invariant counts the number of BPS 
states weighted by a sign. From a more mathematical perspective, the index corresponds 
to topological invariants (e.g. the Euler number or the Betti numbers) of a moduli space of 
objects (of an appropriate category) corresponding to the BPS states. 

One of the seminal papers on BPS invariants of supersymmetric gauge theory on a Kahler 
surface is Ref. [34j by Vafa and Witten. They show that the topologically twisted path 
integral localizes on the instanton solutions, and equals the generating function of the Euler 
numbers of instanton moduli spaces, whose natural compactification is the moduli space 
of semi-stable sheaves. One of their main motivations was to test the strong-weak coupling 
duality [SU] or ^-duality, which acts by SL(2, Z) transformations on the theory. The coupling 
constant g and the theta angle 9 combine to the modular parameter r = J- + ^-duality 



suggests that the generating function of the BPS invariants (3.2) should exhibit modular 



properties if the gauge group is SU(r) or U(r). They tested this in various cases, for example 
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for sheaves with rank r = 1 j 12 j . and r = 2 on P 2 [36| [37J [19]. The generating functions 
for rank 1 were found to be genuine (weakly) holomorphic modular forms. However the 
generating functions for rank 2 transform only approximately as a modular form. These 
functions are (mixed) mock modular forms, i.e. functions which do transform as a modular 
form only after the addition of a non-holomorphic "completion" [41J. 

Ref. [34J has inspired many results in later years. In particular, for r = 2 the dependence of 
the BPS invariant on the polarization J G H 2 (S, Z) was included in the generating functions 
using indefinite theta functions [H]. Moreover, the reduced modular properties for r > 2 
were understood physically as a "holomorphic anomaly" |2"9[ H]. 

Although modularity has proven useful for various computations [291 1391 HI] , physical 
expectations for r > 2 could never be rigorously tested since generating functions for r > 2 
were not known. This was one of the motivations for (57], which computed the generating 
functions of refined BPS invariants for r = 3 on P 2 and its blow-up P 2 , which is the Hirzebruch 
surface Si. A convenient property of Si is that the BPS invariants vanish for certain choices 
of the first Chern class and choice of polarization. Wall-crossing and the blow-up formula 
|38] provide then the invariants in the other chambers and for P 2 J^] 

This article generalizes the computation of the generating function Z r (z, p, r; Si, J) of 
BPS invariants for r < 3 of Ref. [27] to more general Hirzebruch surfaces S^, where —I is 
the self-intersection number of the base curve of S|. The arguments z G C, p G if 2 (S^,C) 
and r G "H in Z r (z, p, r; S^, J) are generating variables for the Betti numbers of the moduli 
spaces, and first & second Chern classes of the sheaves respectively. 



Section 3.1 derives expressions for the generating functions with r = 2 in terms of indefinite 
theta functions [14J and Appell functions of level I [3] [32]. The non-holomorphic but modular 
completed functions Z 2 (z, p, t; S^, J) are determined for zgC (generating function of Betti 
numbers) as well as z = \ (Euler numbers). Due to the presence of these terms the action of 



the heat operator D r on the generating function Z r (p, r; S^, J) (3.1 ) of Euler numbers does 



not vanish, which is known in the physics literature as a "holomorphic anomaly". A novel 



result of the paper is that D 2 Z 2 (p, r; S^, J) in general consists of two terms (3.11 ) 
(1.1) D 2 Z 2 (p, t; S,, J) = C 2 (Imr, J) Z 1 (p, r, S,) 2 + R 2 (p, r; S,, J), 



Refs. |22l I35| computed earlier generating functions for the Euler numbers for rank 3 using different 
techniques. 
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where ^(Imr, J) is a simple function of Imr and J. The appearance of Zi(p, z, r, S^) 2 has 
been conjectured and discussed in the literature before [Ml 1291 0], but the additional term 
i?2(p, -2, t; S^, J) is novel. Remarkably, the additional term vanishes for special choices of J, 
in particular for J = —Kg where Ki is the canonical class of S^. ^\ 

Another important property of the non-holomorphic completion is that it renders Z r (p, z, r; E^, J) 
continuous as a function of the polarization J [25J , which is expected of a physical path inte- 
gral. Although a more intrinsic derivation of the anomaly in physics or algebraic geometry 
is desirable, this gives already important insights. 



Section 3.3 presents the holomorphic generating function for r = 3 (3.12) for r = 3 
and presents the Tables [T]j3] with the Betti numbers for £ = 1. The modular properties of 
Za(p, z, r; S^, J) are much more intricate then for r = 2, and will be discussed elsewhere (5j. 

The outline of the paper is as follows. Section [2] reviews the necessary properties of sheaves 
and Hirzebruch surfaces, including BPS invariants and their wall-crossing. Section [3] de- 
fines the generating functions and gives explicit expressions for r = 1, 2 and 3. The non- 
holomorphic terms and the holomorphic anomaly are determined for r = 2 in Subsection 
3.2[ and for r = 3 Tables with Betti numbers are presented in [373 



2. Sheaves on Hirzebruch surfaces 

The Gieseker-Maruyama moduli space of semi-stable sheaves with rank r on S is the 
natural compactification of the moduli space of instantons with gauge group U(r), i.e. anti- 
self-dual solutions for the field strength: *F = —F. The Chern classes of the sheaf are 
determined by the topological classes of the instanton: 

Most of the following is phrased in the more algebraic language of sheaves, since this setting 
is most suitable for explicit computations. 

2.1. Sheaves and stability. The Chern character of a sheaf F on a surface S is given 
by ch(F) = r(F) + C\{F) + \c\(F) 2 — 02(F) in terms of the rank r(F) and its Chern 



'Note for £ > 2, — K( does not lie in the ample cone of and is therefore not a permissible choice for J. 
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classes Ci(F) and 02(F). The vector T(F) := (r(F), chi(F), cli2(F)) summarizes the topo- 
logical properties of F. Other frequently occuring quantities are the determinant A(F) = 
^)(^(F) - ^^(F) 2 ), and //(F) = Ci(F)/r(F). 

Let C Fx C F 2 C • • • C F s = F be a filtration of the sheaf F. The quotients are denoted 
by Fj = Fj/Fj_! with Tj = T(Fj). With the above notation, the discriminant A(F) is given 
in terms of the topological quantities of Fj and Fj by 

(2.1) A(„ = £ - ^ £ ^™ _ , {Fi) f . 

The notion of a moduli space for sheaves is only well defined after the introduction of a 
stability condition. To this end let C(S) G H 2 (S, Z) be the ample cone of 5*. Given a choice 
J G C(S), a sheaf F is called //-stable if for every subsheaf F' C F, /i(F') • J < //(F) • J, 
and //-semi- stable if fi(F') ■ J < //(F) • J. A wall of marginal stability W is a (codimension 
1) subspace of C(S), such that (/i(F') — //(F)) • J = 0, but (//(F') — //(F)) • J 7^ away from 
M/. 

Let 5* be a Kahler surface, whose intersection pairing on H 2 (S, Z) has signature (1, 62 — 1)- 
Since at a wall, (//2 — //1) ■ J = and J 2 > 0, we have (//2 — //i) 2 < 0. Therefore, the set 
of semi-stable nitrations for F, with Aj > is finite. The ample class J provides natural 
projections c± for an element c G H 2 (S, Z) to the positive and negative definite subspaces 
of H 2 {S, E): 

. , c ■ J J 

(2.2) c + = , c _=c-c+. 

2.2. Some properties of ruled surfaces. A ruled surface is a surface S together with a 
surjective morphism tc : S — > C to a curve C, such that the fibre T, y is isomorphic to P 1 for 
every point y G C. Let / be the fibre of tt, then F2(S,Z) = ZC © Z/, with intersection 
numbers C 2 = -t < 0, f 2 = and C-/ = 1. The canonical class is F s = -2C+(2g-2-£)f. 
The holomorphic Euler characteristic x(^s) is for a ruled surface 1 — g. An ample class is 
parametrized by J mn = m(C + £f) + nf G C(S) with m, n > 0. The following only considers 
surfaces with g = 0, these are known as rationally ruled surfaces or Hirzebruch surfaces. They 
are denoted by and furthermore denotes the canonical class. 

To learn about the set of semi-stable sheaves on for J G C(S), it is useful to first 
consider the restriction of the sheaves on to /. Namely the restriction to E\f is stable if 
and only if F is //-stable for J = J ,i an d in the adjacent chamber [12]. However, since every 
bundle of rank > 2 on P 1 is a sum of line bundles, there are no stable bundles with r > 2 
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on P 1 . Therefore, the BPS invariant Q(T,w; J) (defined in the next subsection) vanishes for 
T = (r(F), -C - af, ch 2 ) with r(F) > 2 and a = 0, 1. 

2.3. Invariants and wall-crossing. The moduli space Aij(T) of semi-stable sheaves (with 
respect to the ample class J) whose rank and Chern classes are determined by T has complex 
dimension: 

(2.3) dime Mj{T) = 2r 2 A-r 2 X {O s ) + l. 

To define the refined BPS invariants f2(r, to; J) in an informal way, let p(X, s) = X^i=o ^ ^i s% 
with hi the Betti numbers 6, = dim if 2 (X, Z), be the Poincare polynomial of a compact com- 
plex manifold X. Then: 

w -A\va c Mj(T) 

(2.4) Q(T,w;J) := p( M j(T),w). 

w — w 1 

The rational refined invariants are defined by [27J: 

tt(T/m,-(-w) m ;J) 



m\V 

See [2H] for a physical motivation of these rational invariants and |2"Uj I3T] for mathematical 
motivations. The numerical BPS invariant fi(r; J) follows from the 0(T, io; J) by: 

(2.5) fi(r; J) = lim (w - u; -1 ) fi(r, u;; J), 

1U— >— 1 

and similarly for the rational invariants Cl(T; J). 

A crucial tool for the computation of the generating functions in Section [3] is the wall- 
crossing formula, which provides the change Afi(r; Jq — > Jc) across walls of marginal sta- 
bility. Ref. [37] gives as criterion for his wall-crossing formula for r = 2 that Kg ■ J < 0, 
which holds for any I and J G C(E^). For r = 3 more complicated wall-crossings appear, 
in particular walls where the slope of three rank 1 sheaves with different C\ become equal. 
Physical arguments suggest that for these walls one could use the wall-crossing formulas of 
Kontsevich-Soibelman (20] or Joyce-Song \17\ since they are shown to hold in both supergrav- 
ity and field theory [S| [2"jlll|. These wall-crossing formulas are derived for Donaldson-Thomas 
invariants, which are defined for 6-dimensional gauge theory on a Calabi-Yau 3-fold [9|. The 
mathematical justification for the use of these wall-crossing formulas for sheaves on surfaces 
is therefore not well established. Ref. [16] gives as criterion for the applicablity that Kg 1 
must be numerically effective (i.e. —Kg ■ D > for any curve in D e H 2 (S, Z)). This would 



exclude the Hirzebruch surfaces with £ > 2. The generating function (3.12) for r = 3 is 
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consistent with the wall-crossing formulas for DT-invariants and in agreement with previous 
results in the literature for £ — 1, but in view of the above requires at least for I > 2 further 
justification. 

Keeping in mind these comments, I continue by giving the explicit change of the invariants 
in case of primitive wall-crossing. To this end, define the following quantities: 

(ri, T 2 ) = rir 2 (/i 2 - Mi) ■ K s , X(Ti, T 2 ; J) = nr 2 (/i 2 - fn) ■ J. 

The change A£7(ri + T 2 , w; Jq — > Jc) for T\ and T 2 primitive is |38l 120] 

(2.6) AQ(T,w;J c ^Jc') = -\ (sgn(X(r 1 , T 2 ; J c ,)) - sgn{Z{T x , T 2 ; J c ))) 

x (w {Fi ' r2) - w - {Ti ' r2) ) n(r u w, j) n(r 2 , w- j). 

with 

sgn(x) = 

The subscript Wc in Jw c refers to a point in C which is sufficiently close to the wall W, such 
that no wall is crossed for the constituent between the wall and Jw c - Note that the wall is 
independent of c 2 . 

For the computation of the invariants of rank 3, one also needs to determine the wall- 
crossing formula across walls of marginal stability for non-primitive charges 21^ + T 2 and 
walls where the slope of three non-parallel charges becomes equal. These can be determined 
using the wall-crossing formulas |20[ ITT] . The result takes a simple form in terms of rational 




invariants and (2.6) |26j. 



3. Generating functions 

This section computes the generating functions of the BPS invariants Q(T,w; J). We 
start by defining the generating functions and a brief discussion of their properties. The 
generating function Z r (p, z, r; S, J) for a Kahler surface S is defined by: 



Z r (p,z,r;S,J) = ^0(I>;J)(-ir ci ^ 



x ^rA(r)- V&p. - i (ci+ri( s /2)^l( Cl +r/f s /2)2 e 2nip-(c 1+ rK s /2) 

with p G H 2 (S,C), w = e 2mz and q = e 2mT . Twisting by a line bundle leads to an isomor- 
phism of moduli spaces. It is therefore sufficient to determine f2(r, w; J) only for c\ mod r, 
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and it moreover implies that Z r (p, z, r; S, J) allows a theta function decomposition: 

(3.1) Z r (p,z,r;S,J) = ^ h r ,„(z,r; S, J) O r , M (p,r; S), 

where the bar over h r ^(z,T; S, J) denotes complex conjugation, and h ryfl (z,r; S, J) and 
< £ > r, f j,(p,T~; S) are defined by: 

(3.2) hr,^r;S,J) = ^ Cl(T, w; J) g rA(r) ~ 12 ^ 1 , 

C2 

QrAp,r;S) = Yl (_l)^ Sg kV2 V ki/2r e 2™p.k 

keH 2 (S,rZ)+rK s /2+fi 

Note that Q r ^(p, r; S) depends on J through k± and does not depend on z. 



The generating function of the numerical invariants £l(T; J) follows simply from Eq. (2.5): 
Z r (p, r; S, J) = lim (w — w' 1 ) Z r (z, p, r; S, J). 



Physical arguments imply that this function transforms as a multivariable Jacobi form of 
weight (|, — |) |34"t 124"] with a non-trivial multiplier system. For rank > 1 this is only correct 
after the addition of a suitable non-holomorphic term [Ml [29] . This is explained for r = 2 



in Subsections 13. II and 13.21 

The functions h r ^ Cl (z,r) and h rfil (r) contain a factor which depends only on the rank r 
and b2{S). It is therefore useful to define 

/r ' Cl(Z ' T) = (^,t)77(t)MS)-i) KcA^r), 



with 9\(z,t) and ^(r) defined by (A.l). The function f r , cl ( T ) follows from / r ci (z, r) by 

Note that the terms of degree < r — 1 in the Taylor expansion with respect to z of f r ,^{z, r) 
vanish. 

A useful relation is the "blow-up formula" which relates the generating function of a surface 
S with that of its blow-up <p '■ S — > S at a non-singular point. Let C\ be the exceptional 
divisor of 0, and take J G C(S), r, and cj such that gcd(r, c\ ■ J) = 1. The generating 
functions /i r>Cl (z, r; 5*, J) and h ryCl (z, r; 5*, J) are then related by [36, EH [38j [23j H3] : 

(3.4) h r ^ cl _ kCl (z,r;S,J) = B t ^{z,t) h rjCl (z,r; S, J), 



8 



JAN MANSCHOT 



with 



B r;k (z,r) 



2^ i 



1 V^ r 2 

2 2^1=1 a i cii—a 



3.1. Rank 1 and 2. This subsection presents explicit expressions for h r , Cl (z, r; S^, J m>n )- 
The result for r = 1 and 5 = is simply [12J: 

/i )C1 (z,t;£^) = 1. 

Note that the dependence on J could be omitted here since all rank 1 sheaves are stable. 
Moreover, there is also no dependence on I. 



n Wj w 2 



o 



m 



FIGURE 1. The ample cone of S 1; together with the three walls for V = 
(2, -C - f, 2), namely for (a, b) = (1, 0), (2, 0), (3, 0). 

To compute the generating functions for r > 2, we use wall-crossing together with the fact 
that Q(T,w, Jo,i) — for c\ = —C + af and r > 2. In the following, cxiE^) is parametrised 
by bC — af. The walls are then at ™ = for r = 2, with m,n > 0. See Figure [lj for the 
walls for A(F) = §, r(F) = 2. One finds [H 123 using Eq. (O): 

f2,c-af(z,r;T l! ,, J m<n ) = -± ^ 0j6eZ |(sgn((26 + l)n - (2a - a)m) - sgn(26 + 1) ) 

( 3> 5) x ^(/-2)(26+l)+2(2o-a) _ w -(*-2)(2&+l)-2(2a-a) J 

x | (2&+l) 2 +i(2&+l)(2a-c*)_ 

These functions are indefinite theta functions j 13 j . which are sums over a subset of the 
positive definite sublattice of an indefinite lattice. Since the sum is only over a subset of the 
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lattice, they transform as a modular form only after addition of a suitable non-holomorphic 
term (depending on f and z) [41] . 

The computation of the invariants for c\ = —af is much more involved since strictly 
semi-stable sheaves do exist for J 0j i or if a = for every J G C(P 2 ). We will circumvent this 
computation by determining the functions f'2,- a f{z, t; S^, J^n) from modular transformations 
of f2,c-af{ z , T 'i ^ii Ji,o)- One can consequently determine the invariants for arbitrary J m n 
by application of the wall-crossing formula. 

We continue by writing f2,c-af{z, t; E^, J m ,n) in terms of two new functions A^t a ^-\{z, r) 
and $™p(z,T): 



with 



f2,pc- a f(z, t; Ht, J m ,n) = A tM {z, t) + $™p(z, r), a,P € {0, 1}. 



€SM = E §(sgn(-(2a-a))-sgn((2&-^-(2a- a )m)) 

a,beZ 

xw (£-2)(2b-f3)+2(2a- a ) g f (2&-/3) 2 +A(2&-/3)(2a-a) _ 



Then Eq. (3.5) gives for A^c-af( z ) T ) for £ > 1 after performing a geometric sumj^j 



e+2 fa(n+l)+n 2(«-2)n 

(3.6) %, 1 ,(^) = r^i T -^ r 



1 . „M(n+l)„,,2(£-2)n 



-.2n+l w 4 



The functions in Eq. (3.6) are specializations of higher level Appell functions [3j |4*2] . 
whose definition is recalled in Appendix A. These functions appeared earlier in mathematical 
physics in the theory of characters of superconformal algebras [10, HHJ [32]. See [33] for a 
recent discussion. This might not be accidental since M = 4 Yang-Mills is well known to 
be related related to 2d conformal field theory by M-theory [21]. Deriving these functions 
explicitly from a 2-dimensional perspective is an interesting direction for future research. 

Analogously to the indefinite theta functions, the Appell functions only transform as a 



modular (or Jacobi) form after addition of a non-holomorphic term. Eq. (A. 3) gives the 



exact expression obtained by [42]. Application of this to our case of interest gives for the 



3 Note that for S^ =0 — P 1 ® P 1 ; the function A ,(o,i)( z j T ) i s undefined while A 0i (i i i)(z, r) = 0. 



10 



JAN MANSCHOT 



completion Atr a p = u(z,T) 



( 



\ 



A 



Z,T) 



A 



l,{a,0) 



Z,T) 



W 1 ,L1 q4.e 



(3.7) 



+ 2-E E 

fc=0 \ ni=2fe + /3< 
\ mod 2 

^2 (sg n (^2) - E {[n 2 + 2(£ + 2)Im z/y) y/y/ij^j w~^ n2 q~^ , 



no — — 2k — a 



with y = Imr and E(x) = 2 J* e~ wu2 du. The four functions Ai^ a ^ transform as a vector- 
valued Jacobi form of weight 1 and index —8 of SX(2,Z) [32| l4"2] . One finds for the action 
of the generators S and T : 

z —1 \ T 



(3.8) 



S : 



T : 



.4 



£,(a,/3) 



r t 



'.z,T + r 



E (-D 



a,/3e{0,l} 



± A 



t,{a,0) 



Z,T 



The modular transformations (|3.8|) together with the single pole in z of the refined invari- 

i fj(r) 3 

i t)(t) 3 



ants (2.4) do fix the functions Ae^ a ^(z,r) to be: 
A W) (z,t) =w 2 2_^ — 



q 2n w A 



9 1 {Az,r) 



A ( 



(0,0) 



q £n' w 2(e-2) 



" + E\ 



q£n 2 w 2(£-2)n 



+ 



g 2 "M7 4 0i(4z,r)" 



This agrees for C\ = f with the generating function in Ref. 



(Corollary 3.4). The 



completion of these functions is given by Eq. (3.7). 



One can show the following relation between A\( a fl\(z, r) and Air a n(z, r) using the quasi- 



periodicity formula (A.4): 
(3.9) 



. . 9 2 (2z,2r) . . . 
A,(i,o)(^r) = ^-^A 1 , (1 , 1) (,,r), 

3 (2z,2t) 

A ' (0 ' 0)( "' T) = ^27) A ' (0 ' 1)( "' r) - 



This relation is understood in algebraic geometry by the blow-up formula (3.4), which relates 
the functions h2, Cl {z,T; Si, J^o) with ci = C — af to those with c\ = —af. For /i2 )Cl (z, r;P 2 ) 



one recovers the result of [4J. The multiplicative relation (3.9) does not hold for £ > 1, since 
S£>! is the blow-up of the weighted projective plane (1, 1,£) at its singular point [8J, and the 
blow-up formula is thus not applicable. 



BPS INVARIANTS OF N = 4 GAUGE THEORY ON HIRZEBRUCH SURFACES 11 

What remains is to complete the indefinite theta functions ^^(z,t). One finds using 
Ref. |H]: 

= EH E ((- 2a + a + 2 ( £ + 2 )Wz/)v / ^) 

a,6eZ 

(3.10) - E (((26 - /3)n - (2a - a)m + 2(2n + + 2)m)Im z/y) sjyjj^j 

xw (e-2)(2b-f3)+2(2a-a) £ (26-/3) 2 + i(26-/3)(2a-a) 

with J^ n = m(£m + 2n). The completion for fific-af follows directly from ji^c-af — 
^,(o,/3) + ^a)?- The non-holomorphic term of the first line in Eq. (3.10) is cancelled by the 



non-holomorphic term of Ai^ a ^)(z, r). Thus for the completion of f2,pc-af (and therefore 
also of h 2 ^c-af) the non-holomorphic part of the second line in Eq. (3.10) suffices. We 
define Z r (p, z, r; S, J) := E M e// 2 (s,,z/rZ) %A Z > r ' 5 ' J ) ^(P> r ! 5 )- 



3.2. Holomorphic anomaly for rank 2. This subsection derives D r Z r (p, r; S^, J) for 
Z) r = <9 T + 4^:^ + and r = 2. Since D r <d r41 (p,r; E^) = for any r, it suffices to determine 
dff rtC1 (r; Ti£, J). For a clear exposition, the generating functions are given in this subsection 
in terms of Kg, J etc. instead of the explicit integers £, m and n etc. 

We determine first the completion /2 iCl (r; E^, J) from the generating functions in the 
previous subsection. The result follows from the following three steps: 



use Eq. ( |3.3 ) after replacing the functions with their completions, 



use that E(z) = 2 /* e'^du = sgn(s)(l - f3i {z 2 )) with z G E and 

poo 

v {x) = / u- u e- nu du, 

J X 

and finally use 



03 (x) = 2x-*e~ nx - 2ir(3i_(x) 



One obtains: 



?2, Cl {r\ £/, J) = /2, Ci (t; S £ , J) + 

CP-C1 \ / 



C 

) 



+ £T 2 (£ £ ,2Z) 
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where c± are given by Eq. (2.2). It is now straightforward to compute 9 f / r)C1 (r; E^, J): 
d f f r;C1 (r;E e ,J) = %K "'J. 3 (-l)^ Ci e 2 ,_ Cl _^(0,r;S,) 

167T V J 2 |/2 



8v^ „~ ~VJ^ 

+ff 2 (E^,2Z) 



After combining this result with 02, C:L (p, t; E^) as in (3.1) and manipulation of the lattice 
sums, one obtains for D 2 Z 2 {p^ 7~; E^, J): 

(3.11) D 2 Z 2 (p,T;T; i ,J) = — %K * J 3 Z x {p,r\ E e , J) 2 



167r y J 2 y 



i 



h lfi (T;Y, e ) 2 Y T c1 (t,I!*)62 jC1 (p,t;I1<), 



vy c 1 eH 2 (T, i ,Z/2Z) 



wherd^] 



T Cl (r,E,) = £ ^.c_^(-l)^g</ 4 g- c2 - /4 . 

* 'J 2 



Interestingly, Eq. (3.11) differs from the conjectured form of the anomaly [34, [291 HI- The 
first line has the expected factorized form, which is attributed to reducible connections or 
polystable sheaves [34J or multiple M5-branes |29j . However, the novel second line does not 
factorize and is less easily interpreted. It does vanish for special values of J, in particular 
for J = —Kf since then Kg ■ c_ = 0. But for i > 2, Kg lies outside C(S) and is thus not a 
permissible choice for J. Viewing the surface as part of a local Calabi-Yau 3-fold geometry, 
J = —Ks corresponds to the attractor point from the point of view of supergravity [25]. It 
is therefore rather interesting that Z(p, r; E^, J) simplifies at this point. 

The function T Cl (r, E^) vanishes also for i = 1 and J = C+f jl], which is not equal to —K\ . 
For this choice, the blow-up formula gives the generating function for P 2 , where J = —K^2 
is satisfied automatically. It is thus in agreement with these examples to conjecture that 
generically for a Kahler surface S, D 2 Z 2 (p, r; S, J) = — Z\(p, r; S, —K s ) if Ks G C(S). 

167T1/2 

Of course, a more intrinsic explanation based on gauge theory or algebraic geometry is 
desirable. 

3.3. Rank 3. This subsection presents the generating functions h^pc-af{z, t; E^, J) with 
(3 7^ mod 3. This condition on (5 ensures that hz } pc-af{ z : T \ E^, J) = for J = J ,i 
analogously to r = 2. The computation of h^c-af{ z i T 'i E^, J) therefore reduces again to 

4 A similar function appeared in Ref. |25| . 
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application of the wall-crossing formula. This is for r = 3 more complicated than for r = 2 
since: 

- the functions hi fix {z, r\ E^, J) do themselves depend on J, and need to be determined 
sufficiently close to the appropriate wall. 

- the total charge T can be of a sum of 3 charges Yli=i Tj such that at a wall W the 
slopes of these three constituents might be equal. This in particular happens for 
"semi-primitive wall-crossing" where T(F) = 2Ti + T 2 . 

Nevertheless, the wall-crossing formulas [20]. [T7] imply a relatively simple form for the gen- 
erating functions I2"T) . One obtains for I > 1: 

f3,(3c-af(z, r; J m>n ) = - ^ §( sgn((36 - 2/3)ra - (3a - 2a)m) - sgn(36 - 2/3) ) 

(3 12) X r w (^-2)(36-2^)+2(3a-2a) _ ly -(^-2)(36-2^)-2(3a-2a)A 

X g ^(3&-2/3) 2 + |(3b-2/3)(3a-2a) 

X f2,bC-af(z, T] T,£, J\3b-2/3\,\3a-2a\), 



for /3 = 1,2 mod 3 and a G Z. Writing out the lattice sums in Eq. (3.12), one finds a novel 
indefinite theta function. It has signature (2, 2) and the condition which determines whether 
or not a lattice point contributes depends quadratically on the lattice vector, whereas pre- 
viously described indefinite theta functions have signature (n, 1) and the condition depends 
linearly on the lattice vector [131 Si] - A detailed discussion of the (mock) modular properties 
of ft,3 )C1 (z, r; E^, J) will appear in a future article |5]. 

Tables [Ip] list Betti numbers for c\ = —C — af with a = 1, 2, 3 and £ — 1, which are 



in agreement with the expected dimension (2.3). One can relate the Betti numbers for 



Ci = — 2C — af to these by using h r ^ Cl = h r _ ci , and /v,ci+k — h r ,a for k G H 2 (S,rZ*). With a 
little more work, one can verify that h^ }C1 {z, r; E 1; J 10 ) satisfies the relations implied by the 



blow-up formula (3.4). 
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c 2 b b 2 b 4 b 6 b 8 b 10 b u bu he h$ hp h 2 h 4 X 

2 1 2 4 4 18 

3 1 3 9 20 37 53 59 305 

4 1 3 10 25 59 119 218 338 450 490 2936 

5 1 3 10 26 64 141 294 562 997 1602 2301 2886 3117 20891 



TABLE 1. The Betti numbers b n (with n < dim c M) and the Euler numbers 
X of the moduli spaces of stable sheaves on with r = 3, C\ = —C, and 
2 < c 2 < 6 for J = (l,e). 



c 2 b b 2 b 4 b 6 h hp b 12 b u b 16 bi S ho b 22 hi h& X_ 

2 11 3 

3 1 3 8 14 17 69 

4 1 3 10 24 53 93 136 152 792 

5 1 3 10 26 63 135 268 470 725 950 1043 6345 

6 1 3 10 26 65 145 310 612 1144 1970 3113 4391 5462 5873 40377 



TABLE 2. The Betti numbers b n (with n < dime M.) and the Euler numbers 
X of the moduli spaces of stable sheaves on Ti 1 with r = 3, Ci = — C — f, and 
2 < c 2 < 6. 



c 2 b b 2 b 4 b 6 b 8 hp bu bu b w h$ b 20 b 22 x 

3 12 3 9 

4 1 3 9 19 31 36 162 

5 1 3 10 25 58 113 192 264 297 1629 

6 1 3 10 26 64 140 288 536 907 1348 1733 1885 11997 



TABLE 3. The Betti numbers b n (with n < dime -M) and the Euler numbers 
X of the moduli spaces of stable sheaves on S x with r = 3, C\ = —C — 2f, and 
3 < c 2 < 6. 



BPS INVARIANTS OF N = 4 GAUGE THEORY ON HIRZEBRUCH SURFACES 15 

Appendix A. Modular functions 

Define q := e 2mT , w := e 2mz , with r G HI and z G C. The Dedekind eta and Jacobi theta 
functions are defined by: 



71=1 

„2 



(A.l) 0i(s,r):=i (-irk^w r , 



i 



2 

,2 



n /2 w n . 



9 3 (z,r) :=J2« r ' 
The Appell function at level £ is defined by: 

(A,) ^,., T) = ^ E (=!g^ , 

with a = e 2 ™ and 6 = e 2 ™. In order to give the completion Aj>(u, v, r), define 
R(u,t) = y sgn(r) — E (jr + lmu/y)y/2y^ 



x(-l) r -3 a -V r2/2 , 
with = 2 e~™' ' du. The completion Ai(u,v,t) is then given by 

. £-1 

(A.3) r) = A e {u,v,r)+ l -J2 a k O^v + kr + {£ - l)/2,£r) 

k=0 

x - v -kr — (£- l)/2, £r), 

/ -£ 1 \ 

and transforms as a multivariable Jacobi form of weight 1 and index | ( ^ ) ' ^' ie 
pell function for i = 1 is related to the Lerch- Appell function: fi(u, v, r) = Ai(u, v, t)/$i(v), 
which satisfies the quasi-periodicity property [41J: 

(A.4) fi(u + z, v + z, t) -fi(u,v,r) = ^tL^ ^ ~j~ ~ 4^T~ — V' 

r) 0(t>, r) 0(it + 2, r) 0(v + z, r) 

for u, v , w + z, v + z ^ Zr + Z. 
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